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Abstract 

This short note presents the derivation of a new a priori estimate for the Oldroyd- 
B model. Such an estimate may provide useful information when investigating the 
long-time behaviour of macro-macro models, and the stability of numerical schemes. 
We show how this estimate can be used as a guideline to derive new estimates for 
other macroscopic models, like the FENE-P model. 



1 Introduction 

We consider the Oldroyd-B model: 

I du \ 

Re ( — +u.Wu) = (l-e)Au- Vp + divx, (1) 

div (u) = 0, (2) 

+ u.Vt = Vut + t(Vu) t - — r + -^(Vu+ (Vu) T ) , (3) 
at We We 

where the Reynolds number Re > 0, the Weissenberg number We > and e € (0,1) 
are some non-dimensional numbers. We suppose that the space variable x lives in 
a bounded domain V of R d . This system is supplied with initial conditions on the 
velocity u and on the stress tensor r. For simplicity, we assume no-slip boundary 
conditions on the velocity it: 

u = on dV. (4) 

We suppose that the initial data and the geometry are such that there exists a unique 
regular solution to (J])-© and our aim is to derive some a priori estimates on this 
solution. 

Let us introduce the so-called conformation tensor A = + Id. The partial 

differential equation (PDE) on r translates into the following PDE on A: 

8A 1 1 

— +u.\7A = \7uA + A{X7uf-— A + — Id. (5) 
at We We 

One can check that if 

We 

A(t = 0) = r(t = 0) + Id is a positive definite symmetric matrix, (6) 

then this property is propagated forward in time by (and, in particular, r is 
symmetric). Assuming uniqueness of solution, this can be proven for example by 
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using the probabilistic interpretation of A as a covariance matrix, as explained in 
Section [31 We will assume throughout this note that is satisfied. Concerning the 
importance of positive-definiteness of A, we refer for example to Section 9.8.10] 
and also to the recent work [3j[4]. 

In Section [21 we recall how the classical a priori estimate for the Oldroyd-B model 
is derived. Next we show how it can be used to derive some bounds on the stress 
tensor, provided the initial condition satisfies detA(f = 0) > 1. In Section [3j we 
establish a new estimate, which comes from an entropy estimate on the micro-macro 
model associated with the Oldroyd-B model (see |5]). This estimate provides bounds 
on (u,t) without any assumption on r(t = 0) (apart from ([6])). This new estimate 
could be useful to study the longtime behaviour of some macro-macro models, or to 
analyze the stability of some numerical schemes. Current research is directed towards 
clarifying this. 



2 The classical estimate 

Let us first introduce the kinetic energy: 



2 

We easily obtain from (J])-((2]): 



E{t) = \ I |«| 2 . (7) 



v 



Re-^=-(l-e) / |Vw| 2 - f r : Vtt 
dt Jx> Jv 



(8) 



where for two matrices A and B, we denote A : B — A^jBij = tr (A T B). On the 
other hand, taking the trace of the PDE |(3]) on t and integrating over V, we get: 

d P n P „ 1 , 

trr = 2 / Vu:t- — / trr. 



dt J v J v We Jv 



We thus obtain the following estimate: 



s(fX'«" 2+ U trT 



(!>) 



Remark 1 In terms of A, the energy estimate (0) writes: 
d ( Re f , , 2 £ 



, \u\ z + / txA 

dt V 2 ,/p 2 We J v 



f |Vm| 2 + ^^ { tr (A - Id) =0. 
'Jv 2We 2 Jv 



(10) 



In Lemma Q] below, we prove that trr is positive if det A(t = 0) > 1. This result 
combined with the estimate J9J thus yields some a priori bounds on (it, r) provided 
det(A)(t = 0) > 1. In particular, it shows that u and r go exponentially fast to in 
the long time limit, using (J9j) and the Poincare inequality: J v \u\ 2 < C |Vtt| 2 . 

Lemma 1 Let us assume that det A(t — 0) > 1. Then, we have Vt > 0, det A(t) > 1 
and this implies that tr t{£) > 0. 

Proof: Using ((5j) and the Jacobi identity (which states that for any invertible matrix 
M depending smoothly on a parameter t, ^ In det M = tr (-^ _1 ^|f)), we have: 

91n( f A) + «.V ln(dct A) = -Ltr (A^ Id) . (11) 
at We 
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Since for any symmetric positive matrix M of size d x d, 

(detM) 1 ^ < (l/d)trM, (12) 

we obtain 

^MdetA) + u v ^_ / /d _ x 

OT We V / 

which we can rewrite in terms of y = (det A) 1 /'': 

We f|| + u.Vy) > (1 - y) . (13) 

This shows that y > 1 if y(t = 0) > 1, and thus that det A > 1 if det A(t = 0) > 1. 

Indeed, using the characteristic method (by integrating the vector field u(t,x)), 
one can rewrite l|13p as 

We^| > (1 - y) . 

Now, if y does not remain greater than 1, consider the first time t Q such that yito) = 1. 
We have on the one hand ^f(to) < and, on the other hand (1 — y(<o)) = 0. We 
reach a contradiction. 

We thus have det A > 1 and therefore, using again lf!2|) . trA > d. Since r = 
^g(A — Id), this is equivalent to trr > 0. 

Remark 2 7/ det A(t = 0) < 1 (which is the case iftrr(t = 0) <Q), Equation figj) 
shows that det A grows along the characteristics as long as det A < 1 . 

3 Entropy estimate 

We now consider a micro-macro (or multiscale) formulation of the Oldroyd-B model 
and some estimates based on entropy, inspired from [5]. 

3.1 General derivation of the entropy estimate for micro-macro 
models 

We consider the following system: 

( du \ 
Re -7T-(t, x) + u(t, x).Wu(t, x) = (1 - e)Au(t, x) - Vp(t, x) + div r(t, x), 



dt 

div (u(t,x)) = 0, 

r(t, x) = 4r\ I ( X ® Vn(JT))^(t, a;, X) - Id) , 
We VJr^ / 

^(t,sB,jr) + ti(t,»).v Bl v(*,aj,-x:) 

at 

= -div x ( (v x u(t, x)X - VII(X)) *)) + A *V'(*, x, X). 

(14) 

This system is supplied with initial conditions on the velocity u and on the distribu- 
tion ip. We recall that we suppose no-slip boundary conditions (j4]) on the velocity u. 
This system corresponds to a micro-macro model of polymeric fluids, the polymer 
being modelled by two beads linked by a spring with potential energy II. The con- 
figurational variable X e R d models the end-to-end vector of the polymer. For more 
details on the modelling, we refer to [DIB]. 

Notice that we could rewrite the former system as a system coupling a PDE and 
a stochastic differential equation (SDE), replacing the last two equations by: 

r(t, x) = ±- (E (X t (x) <g> WU(X t (x))) - Id) , (15) 
dX t (x) + u(t, x).V x X t {x) dt 

V x u(t, x)X t {x) - JL-Vn.(X t (x)) \ dt + -^=dW t . (16) 
2We / VWe 
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There, E denotes the expectation, W t denotes a d-dimensional standard Brownian 
motion independent from the initial condition (Xq{x)) x( zt> which is such that, Vcc € V, 
the law of X (x) is ip(0, x, X) dX. 
Let us introduce the kinetic energy: 



m =\L 



We easily obtain: 
dE 



Re 



dt 



-(1-s) 



|V«| 



e 

We 



/ / (X® Vn(X)) : VwV- 

J I? jR d 



We now introduce the entropy of the system, namely: 

fi>(t, x ,xy 



H(t)= f f ^(t,x,X)ln 

= 11 U(X)^j(t,x,X)+ f f ^(t,x,X)\n^(t,x,X))+C\ 



(17) 



(18) 



(19) 



with 



expj-npQ) 
J R ,exp(-n(X)) : 



(20) 



and C = rn(/ Rd exp(— II(X)))|2?|. The function H is actually the relative entropy of 
tp with respect to the equilibrium distribution ■f/'oo- 
After some computations (see [5]), we obtain: 



dH 
dt 



2We 



V JR d 



Vln 



/ / (X<g>VII(X)) : Vm?/ 1 - (21) 



Therefore, introducing the free energy F(t) = E(t) + ^^H(t) of the system, we 
have: 



(22) 



Using a logarithmic Sobolev inequality with respect to V'oo and a Poincare in- 
equality for u G Hq(D), one can then obtain exponential convergence to equilibrium 
lim t ^ 00 (it, ip) = (0, ^oo) (see [5]). For some generalizations to the case u / on dT>, 
we refer to |5|. 




3.2 The Oldroyd-B case 

Let us consider the Hookean dumbbell model, for which the potential II of the entropic 
force is: 

|| x\\ 2 

U(X) = U-JL. (23) 

By Itd's calculus, it is easy to derive from ifTSj) that A = E(X t ® X t ) satisfies the 
following PDE: 

8A 1 1 

— +u.\7A = \7uA + A{X7uf-— A + — Id. (24) 
at We We 

This translates into the following PDE for r = ^^(A — Id): 

^ + u.Vt = Vut + t(Vu) t - -^-t + (Vu + (Vtt) T ) . (25) 
at We We 
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The Hookean dumbbell model is thus equivalent to the Oldroyd-B model (at least for 
regular enough solutions). 

If x, .) is Gaussian (with zero mean), so is ip(t, x, .): 



iP(t,x,X) 



X T A~ 1 X 



(2n) d / 2 ^/M{A 



: exp 



where A = E(X t (g> X t ) = J R d X ® X ip(t, x, X) dX denotes as above the covari- 
ance matrix of X t , which depends on time and also on the space variable x. The 
covariance matrix A is symmetric and nonnegative. Moreover, since for almost all 
t > 0, J v j Rd ip(t, x, X) In ( ^''^x^ ) < oo, then for almost alH > and for almost 
all x e T>, A is positive. 

The following explicit expression of the relative entropy can then be derived: 



/ / ip(t,x,X)]n 

Jv JR d 



f ^(t,x,X) 

V ^oo(X) 



dX = [ -(-ln(detA)-d + trA). 
Jv 2 



On the other hand, 



/ / i>{t,x,x) 

JV JR d 



V x ln 



( j>{t,x,X) 
V MX) 



dX 



tr ((Id - A 



-1\2 



Rewriting (|22|) . we thus obtain the following estimate, in terms of A: 




(26) 



This is, in the specific case of Hookean dumbbells (that is Oldroyd-B model) the 
macroscopic version of (|22|) . 

Since — ln(det(A)) — d + tr (A) > 0, this energy estimate yields some a priori 
bounds on (it, A), and thus on (it, t). In sharp contrast to the classical estimate 
it provides bounds on (u,t) without any assumption on r(t = 0) (apart from ©). 
Using a Poincare inequality and the factlll that, for any symmetric positive matrix M 
of size d x d, 

-ln(detM) -d + tvM < tr ((Id - M^fM) 

exponential convergence to equilibrium (lim t _ >00 (tt, A) — (0, Id)) can be obtained 
from 11261). 



Remark 3 Notice that 



can be schematically obtained as fitfl 



2 We 



V 



Remark 4 lfip(0,x, .) is not Gaussian, it is always possible to replace it by a Gaus- 
sian initial condition with the same mean and variance, so that the macroscopic quan- 
tities (u,p,A) would be the same for the two initial conditions. 



3.3 Application to related macroscopic models 

The energy estimate ([26| can be used as a guideline to derive energy estimates for 
other macroscopic models, even though they cannot be recast as a microscopic model 
of the form (14]) . 

1 which can be seen as the logarithmic Sobolev inequality for Gaussian random variables translated on 
their covariance matrices 
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Let us consider the example of the FENE-P model O E] , for which 



We \l-tt(A)/b 



dA 



— + u.VA = VuA + ACVuY 
at 



-Id 



1 



Wel-tr(A)/6 We 



(27) 
(28) 



For this model, we assume J6|), and also that tr (A)(t = 0) < b, and this property is 
propagated forward in time by (|2"8| (see [6]). 

Using the same ideas as for the Oldroyd-B model, we consider the "entropy" H(t) — 
— ln(det A) — i>m(l — tr (A)/b), and we compute its time-derivative: 



d_ 

dt 



-61n(l - tr (A)/b) 



Vm : A 



1 



tr(A) 



d_ 

dt 



p l-tr(A)/6 WeJ v \ (l-tr(A)/fo) 2 



l-tr(A)/&y ' 
(29) 



/ ln(det(A)) = — [ ( ^rr t-7 +tr(A" 

7p 1 1 " WeJ v \ l-tr(A)/6 1 



-tr 

Combining these expressions with JH) , we obtain 



(30) 



dt 



Re 



2We 



|V«| 



/ (-ln(detA)-61n(l-tr(A)/6)) 
tr(A) 



2r/ 



2We 2 Jv V(l-tr(A)/&) 2 l-tr(A)/6 



tr fA" 



0. 



One can check that for any symmetric positive matrix M of size d x d: 

b 



and that 



ln(det(M)) - b In (1 - tr (M)/6) > -(6 + d) In 



ln(det(M)) -61n(l -tr (M)/b) + (b + d) In 



b + d 



> d 



b + d 



< 



tr (M) 



2d 



(l-tr(M)/6) 2 l-tr(M)/6 



+ tr(M _1 ) 



(31) 
(32) 

(33) 
(34) 



The proof of these inequalities is tedious and can be done by diagonalizing the ma- 
trix M. 

Equation (|32|) shows that 

t I |u|2 + 24 L (" ln(det A) ~ hln (1 - tr + (5 + d) 111 (rh) 



is a non-negative quantity, and thus that (|31|l indeed yields some a priori bounds 
on (u, A). 

Equation l!34tl (which plays the role of the log-Sobolev inequality in the micro- 
macro models) shows that the estimate lf3Tjl can be used to prove exponential conver- 
gence to equilibrium. 
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